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Abstract 

Inclusive single jet production in hadron collisions is considered. It is shown that the QCD 
parton model predicts a nonmonotonic dependence of the inclusive cross section on the 
fraction of the energy deposited in the jet registered, if it is normalized on the same cross 
section measured at another collision energy. Specifically, if the cross section is normalized 
by the one measured at a higher collision energy, it possesses a minimum which depends on 
jet rapidity. This prediction can be tested at the Fermilab Tevatron, at the CERN LHC, 
and at the Very Large Hadron Collider under discussion. 
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The parton model, improved with the QCD running of the coupling and of the structure 
functions, has the well-deserved status of a paradigm. Its diverse justifications are of two 
kinds: the qualitative justifications of the early days of the parton model (like the prediction 
of the relation between structure functions Fi and F 2 governed by the spin of the partons 
(Gross relation), for a review see, e.g., Ref. |L]) and the present-day abundant quantitative 
justifications (for a recent review, including a list of structure function parameterizations, 
see Ref. §). 

The advent of the Fermilab Tevatron and the CERN LHC provides a new testing ground 
for the parton model — the kinematic conditions when the energies of the produced hadrons 
are large enough to be described by perturbation theory and, at the same time, are much 
smaller than the total energy of the collision (semi-hard kinematics). Because the parton 
model was originally invented and subsequently tested for the hard kinematics, the second 
condition makes it plausible that a substantial modification of the parton model will be 
needed to describe this semi-hard kinematic region. Note that for hard kinematics the 
transverse energy of the produced hadrons is on the order of the total collision energy. In 
particular, some fusion of Regge theory (see, e.g., Ref. ||) and perturbative QCD may be 
required for the semi-hard kinematics. The Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
Q and, more generally, the experience with the resummation of the large energy logarithms 
(for a review see Ref. 0) may be useful in constructing this new model. 

The range of applicability of the QCD improved parton model is a subject of controversy 
at the moment. There are statements (see, e.g., |], [71) that the fitting capacity of the con- 
ventional parton model is sufficient to accomodate all the data on parton structure functions 
available at the semi-hard kinematics. On the other hand, the data from HERA on forward 
jet production at small-x || may be interpreted as a manifestation of the BFKL Pomeron 
0, which seems to be beyond the limits of the conventional parton model. The situation is 
further complicated by the observation that the range of applicability of the parton model 
may be different for different observables. In particular, the cross sections of processes with 
specific kinematics exhibit breakdown of the applicability of finite order perturbative QCD 
via the loss of insensitivity to the choice of the normalization scale. On the other hand, some 
dedicated combinations (ratios) of cross sections may be less sensitive to the inclusion of the 

since, 



higher order corrections. An example is the scaled cross section ratio flO| , 11, 12, 12 
as it follows from Ref. [113, it is relatively insensitive to the inclusion of the next-to- leading 
(NLO) correction. 

Under these circumstances, it is crucial to have qualitative predictions from the conven- 
tional QCD-improved parton model (without resummation of the energy logarithms) for the 
new kinematic domain. If the predictions would turn out qualitatively incorrect, a substitute 
for the parton model would become indispensable. 

In this paper, we present such a prediction. It is a prediction for the ratio of inclusive 
single jet production at a smaller energy a/J/v of the hadron collision to the one at a higher 
energy JJ^: 



R ( x y \ _ f s N d(x(s N ) \ /f s D da(s D ) \ 
V dxdy )' V dxdy J 

Here the cross section is made dimensionless by the rescaling with the corresponding total 
invariant energy of the collision squared sn(sd)- The ratio depends on the (pseudo)rapidity 
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y = l/21n(fc-|_/fc_), where k± = E ± k$ are the light-cone components of the momentum of 
the produced jet, and on the fraction of the energy x = (k + + k_)/ -y/Si, i — N,D deposited in 
the jet produced (sjv is used for the definition of x in the numerator, in the denominator, 
so x varies from zero to unity for both energies) . Note that this scaling variable coincides in 
the center of mass system with xr = E/E max = 2E/^/s , the radial Feynman variable. 
The ratio R (taken at y ~ 0, i.e., for jets perpendicular to the collision axes) was used 



in Refs. (JTOJ, [11], |T2| as a means to test QCD predictions for scaling violations. Note that 
without scaling violations the ratio R is exactly unity. The x-dependence of R comes from 
the presence of Aqcd m the running coupling and in the parton distribution functions. 



Perturbative QCD calculations of Ref. with hard kinematics (Q 2 ~ s) predict for R 
at y = a steep increase around the value of 1.8-1.9 for x growing in the range above 0.1 (for 
the case y / Siv/y / so = 0.63 TeV/1.8 TeV ). For moderate x, the prediction is in reasonable 
agreement with CDF data 111], |I2|. For x < 0.1 calculations are above the preliminary data 



of CDF |[L2|| . This was among the reasons for the conclusion of Refs. |14| , p~5f| that NLO 
perturbative QCD |L6| with hard kinematics is insufficient for the description of absolute 
cross section of jets with transverse energy less than 50 GeV within accuracy 10%. It was 



shown in Ref. [HJ that resummation of the energy logarithms restores the agreement between 



theory and experiment. 

In this paper we present the following result: the QCD-improved parton model predicts 
that R is not a monotonic function of its arguments, i.e. the single jet production cross 
section, if measured in the natural units of the same cross section taken at another (higher) 
energy of the collision, has extrema. Namely, it has minima ("dips"): there is a value of x for 
each y with the smallest ratio of jets produced. The reason this fact was overlooked is that 
for y = (the only value for which the calculations were reported earlier) the minimum is at 
a value of x too small to be inside the acceptance of the existing detectors (xdi P (y = 0) < 0.01 
at the Tevatron). 

Fig. 1 presents the ratio for energies 0.63 TeV/1.8 TeV at the Fermilab Tevatron, Fig. 
2 for energies 6 TeV/14 TeV at the CERN LHC, and Fig. 3 for energies 6 TeV/100 TeV 
of the Very Large Hadron Collider (VLHC) proposal fi~8f . Each curve on the plots presents 



the dependence of the ratio on x at different values of rapidity y. Each curve ends at a 
lowest value of x where otsiQ 2 ) — 47r/[(ll — 2/3nj) \n(Q 2 / Aq CD )} has the value of about 
0.5 (it corresponds to Q = 0.7 GeV, and Q was taken to be half of the transverse energy of 
the jet produced). For lower values of x perturbative theory becomes unreliable because the 
coupling approaches unity. 

There is another lower bound on the values of x at which our plots make sense, because 
there is a lowest energy for which the jet may be resolved. This energy is accepted now to be 
around 5 GeV ^ which corresponds to x > 0.016 for the Tevatron (Fig. 1), and to x > 0.0017 
for the LHC (Fig. 2). 

The cross sections for the plots of Figs. 1-3 were calculated with the formula 



dydx x 



4 At the HERA, lepton-hadron collider, jets are resolved from E± = 3 GeV, and at the Tevatron, hadron- 
hadron collider, e.g., CDF Collaboration used to tagging jets from Ej_ = 8 GeV p9|. 
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where Ca = 3 for 577(3), and the two terms in the brackets of the right-hand-side correspond 
to the contributions with balancing jet of rapidity over (I(x, y, Q 2 )) and below (I(x, —y, Q 2 )) 
the one of the registered jet (the balancing jet is the unregistered jet whose transverse 
momentum balances the one of the registered jet; the presence of two terms corresponds 
to the symmetry of the cross section by which it is an even function of y); I(x,y,Q 2 ) is a 
convolution of the parton distribution functions; the normalization point for the coupling and 
for the parton distribution functions is chosen to be Q = 0.5E±, where E± is the transverse 
energy of the jet produced 0. 




In the region of small x (the only one we are interested in), the convolution I(x,y,Q 2 ) 
of the parton distribution functions can be expressed (if one neglects powers of small x) via 
the effective parton distribution functions F(x, Q 2 ) and the effective subprocess of Ref. [pOfl : 

I(x, y, Q 2 ) = / dzF(<P + (x, y; z), Q 2 )F(0_(a;, y; z), Q 2 ), (3) 

where 

M^y,z)= { + eT2y , ( 4 ) 

and z max , z min are determined by 4>+(x,y; z max ) = 1, 4>-(x,y; z min ) = 1. In the calculation, 
F(x, Q 2 ) of the proton was used which coincides with that of the antiproton. Thus, the plots 
of Figs. 1-3 correspond to pp or pp collision. 

There is an important issue concerning the accuracy of present leading order (LO) cal- 
culation. The most important advantage of the scaled cross section ratio is that this is the 
ratio of two perturbative series with the same coefficients and with different scales in the 
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Figure 3: "Dips" at VLHC energies 



running coupling. These scales are defined by the two initial collision energies at fixed scal- 
ing variable. One can show that the theoretical accuracy of the ratio in LO of perturbative 
QCD is at least not less than the accuracy of NLO calculations for absolute cross sections. 

The minima in Figs. 1-3 originate from a competition between the running of the parton 
distribution functions and the running of the coupling constant. Namely, the ratio with 
frozen parton distribution functions is decreasing monotonously (this tendency is realized at 
small x), while the one with frozen coupling constant is growing monotonously (this latter 
tendency is realized for x larger than the position of the minimum). 

We suggest the following potential implications of the minima we have predicted with 
the parton model: (i) If one observes the minima experimentally, one employs the orthodox 
parton model and tries to account for observed positions and depths of the minima by taking 
into account higher order corrections, in particular, resummation of the energy logarithms, 
(ii) If one does not observe the minima experimentally, more radical changes are motivated 
such as an alternative model of the elementary constituents inside the hadrons for semi-hard 



asymptotics. One example might be the color dipole model [21 



Finally, we comment on the possibility for searching the minima at the Fermilab Tevatron, 
at the CERN LHC, and at the VLHC: Positions of the minima for the Tevatron energies 
(see Fig. 1) seem to be reached for both D0 and CDF detectors. The minima of the LHC 



plot (see Fig. 2) seem to be well inside the acceptance of, e.g., the FELIX E3|, the ALICE 



and the CMS detectors. 
We take the ratio of 6 TeV/14 TeV for the LHC, because, in addition to 14 TeV pp 
collisions, lead-lead collisions at the LHC are planned with the collision energy of 6 TeV per 
nucleon-nucleon collision. Since nuclear collisions bring in nuclear effects which can distort 
our predicted curves, we also considered the ratio 6 TeV/100 TeV (Fig. 3). These latter 



predictions correspond to the VLHC energies [fL£ 



Further consideration should be given for which pair of energies and value of rapidity 
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are most convenient for an experimental search of the minima. Also, more work is needed 
to make quantitative predictions for the locations and the shapes of the dips with the NLO 
corrections taken into account. 

It is worth noting that in the case of nuclear collisions the effects of initial nuclear parton 
distributions (small-x EMC-effect ||25|| ) and dynamical effects like quark-gluon plasma, jet 
quenching, etc. |26j will demand special consideration. We plan to consider these and related 
issues in future efforts. 

Before making conclusions, we would like to note that many of the above ideas can be 
studied also for the case of heavy quarkonium production, where similar phenomena should 
be present [F27 . 



To sum up, we find a new qualitative prediction of the QCD-improved parton model for 
hadron collisions and suggest its use to test the applicability of the parton model for certain 
regions of high energy hadron collisions. 
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